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ABSTRACT 



A method has been devised for integrating the pressure-scaled water- 
vapor and CO 2 absorber masses from the surface to the O - levels used 
for prediction in the Naval Postgraduate School primitive equation model. 
Using empirical expressions for atmospheric absorptivity, scattering- 
reflectivity, cloud-reflectivity and earth-surface reflectivity, the 
useful solar insolations absorbed at earth and in the key atmospheric 
layers above the earth have been formulated. The terrestrial cooling 
effect at earth and in these same key atmospheric layers have also been 
formulated using recently published emissivities for the joint effects 
of water vapor and CC^. As a result, the radiative heating (cooling) at 
key levels may be ascertained for input into the atmospheric thermodyna- 
mic equation. In addition, application of atmospheric boundary-layer 
modelling permits determination of the surface-layer turbulent transports 
(in the vertical) of sensible and latent heat at the earth 1 s surface. 
These turbulent transports act in such a way as to provide a heat balance 
at the earth's surface. 
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1. Introduction 



The radiation package described here is primarily an outgrowth of that used 
in the UCLA general circulation model (Gates et al. 1971), although some inno- 
vations from the corresponding model of NCAR (Oliger et al. 1970) and of the 
Geophysical Fluid Dynamical Laboratory (Manabe & Strickler 1964) have also 
been introduced. 

The description here of the radiation package and its bearing upon the heat- 
budget of the earth and atmosphere is similar to that in use by Fleet Numerical 
Weather Central (Kesel & Winninghoff 1972) in application to the FNWC primitive- 
equation (P.E.) forecast model. The most obvious similarities are in the identi- 
cal use of the same five CJ-levels (with Q = p/TT) for purposes of data-input, and 
also of similar criteria for the occurrence of a one-layer cloud deck. The main 
difference in the present from that of the FNWC model lies in the inclusion of 
water-vapor and CO^ reduced optical masses in the o-layer 0.4 to 0.2. This in- 
clusion is made possible by a realistic method of extrapolation of the water- 
vapor specific humidity to levels above CJ=0.4(see Appendix A) . A more clearly 
defined system of integrations for the computation of net flux of terrestrial 
radiation is employed in this study (compared to that of FNWC) by means of new 
empirical absorptivities of water vapor and CO^ (Sasamori, 1968). 

The vertical-structure parameters needed to specify the properties of the 
radiating-absorbing model atmosphere sounding are given in Section 2. The 
solar-radiation package is described in terms of the model in Section 3, and 
the terrestrial radiation package is developed in Section 4. 

2. The radiative vertical-structure parameters. 

(a) Computation of absorber masses. 

The five-layer model of Fig. 1 is similar to that in use in the FNWC 
primitive equation model. Since the half-layer O - levels, denoted by the 
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Level - 
identilers 



O- level 



Layer (0,2) 



(2,4) 



(4,6) 



( 6 , 8 ) 



( 8 , 10 ) 



0 0=0 

1 V,T,q 0.1 

2 a=0.2 

3 VTq 03 

4 0=0.4 

5 VTq ---0.5 

6 0 = 0.6 

7 VTq 07 

8 0 = 0.8 

9 V,T,q 0.9 

10 0 = 1.0 



Fig. 1. Definition of the water vapor and CC^ absorber mass parameters in 
the model. Moisture and temperature parameters are introduced at the odd- 
levels 10. Oo as in the analysis scheme. Levels are identified by their 
values on the 10.0a scale, while layers are identified by their level 
boundary indices [e.g. (8,10)], in parentheses. 
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dashed lines are the primary levels for data-analysis of W, T, and q 
(= mixing ratio), there are a total of 10 levels important for radia- 
tional computations in the model. This happens because for each of 
the five levels where q is known, it is possible to determine the 
pressure-scaled water vapor absorbed mass AU(2k-l) and (similarly that 
of AC (2k-l ) ) , associated with the center of each of the indicated 

five odd numbered levels. The appropriate formulas are 



AU(2k-l) = 



l2k-l Ap ( p 2k-l r 72 



(2-1) 



for water-vapor absorber mass in the layer 2k-2 to 2k. Here Ap = . 2TT 

is a constant of the analysis scheme, and q^^. ^ is the mean value of the 

water vapor mixing ratio centered in the layer (2k-2, 2k). The pressure- 

0.72 

ratio factor (p /p ) is the proposed scaling factor for collisional 

ZK— .L O 

line-broadening, after Moller (1964). This formulation due to Moller, 
disposes of the necessity of temperature scaling in the usual Lorentz 
line-broadening formula (c.f., Danard, 1969), where the pressure-ratio 

scaling is linear and a temperature-ratio scaling is accomplished by 

T - . 5 

(— ) # times the mass within a layer, 
o 



For Moller also proposes pressure-scaling similar in form to (2-1) 

but with a pressure-ratio exponent of 0.65. For simplicity a single 
pressure exponent of .72 will be adopted for both water-vapor and CO^ 
absorber masses. Thus the layer-reduced absorber mass of CO^ in N.T.P 
cm per cm^ column is 

-4 Ap ^2k-l 

AC(2k-l) = 3.14 x 10 4 ^ ( ■ ) (2-2) 

SPo p o 
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In both Eqs. 2-1, 2-2 the subscript "o" denotes an N.T.P condition, e.g., 

= 1013.25 mb, and is the standard density at pressure p = P q and 

T = 273 . 16°K. 
o 

For computation of radiative transfer, the total absorber mass (pro- 
perly pressure-scaled) above a designated pressure level is needed. For 
our purpose here, it is most convenient to sum all absorber masses rela- 
tive to the earth* s surface, where 2k = 10. Thus the reduced absorber 
masses of Eqs. 2-1, 2-2 when applied to the layer (8, 10) become 



^9 ^9 72 

“< 9 > = u 8,io-r Ap ( r>' 

O 

Ac (9) = c = 3,14 x 10 ~ 4a P ( li )- 72 

* C(9) C 8,10 go o ( p Q ) 



(2-3) 



The dual notation of (2-3) suggests the further procedure to be followed 
for incrementing absorber mass over specified layers in the atmosphere. 
The notation of (2-4) indicates the procedure for incrementation of mass 
adopted 



6,10 
' „ 


" U 8, 10 + 

= C - f - 


U 6,8 
C „ 


6,10 


8,10 


6,8 



(2-4) 



Here the method involves incrementing absorber mass only upward over 
successive layers in the five-layer atmosphere. The full atmospheric depth 
of absorber mass is likewise denoted by the symbolism 



U 0,10 " £ =5 u 2k-2,2k 



o 

C 0,10 = ^ =5 C 2k-2,2k 



(2-5) 
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where each term 2k is similar in form to U g with the subscript 

9 replaced by 2k-l. 

Since it is understood that the integration of absorber-mass proceeds 
only from the surface upward, the second subscript "10" in (2-5) will be 
deletedo In line with these concepts, the following table summarizes the 
parameters (and their elevations) needed at each gridpoint to perform a 
radiation calculation. 

Table 1. Notation for the distribution of radiative 
properties at even levels over gridpoints in the five- 
layer model. 



a = p/tt 


1.0 


0.8 


0.6 


0.4 


0.2 


0.0 


10a 


10. 


8. 


6. 


4. 


2. 


0. 


T a 


T 10^s 


T 8 


T 6 


T 4 


T 2 


T o 


U a 


u uf 0 


U 8 


U 6 


U 4 


U 2 


u o 


c a 


C 10- 0 


C 8 


C 6 


C 4 


C 2 


c o 



Since the temperature and humidity analyses are carried along only at 
the odd levels a =.9, .7, .5, .3, .1, it is generally necessary to extra- 
polate a temperature at the surface, using 

T 10 « T s = -5(3T 9 - T ? ) (2-6) 

The temperatures at the other even-subscripted levels are available by 
linear interpolation between temperatures at odd-levels of 10a. In general, 
no effective temperature at level "0" is available, and in fact the topmost 
available temperature analysis gives T^ , which will in general be taken 

‘''This will continue to be the case unless stratospheric temperature- and 
moisture-modelling are undertaken in future. Stratospheric modelling 
should be highly correlated with tropospheric features below. 
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as the effective temperature of the topmost absorber-mass increments of 
water vapor and CO ^ (Uq and Cq), respectively, 

2(b) Vertical extrapolation of the water-vapor mass to upper levels. 

Smith (1966) has shown the existence of a strong correlation between 
q(p) and the surface value through the extrapolative formula 

aM = ( £_ / ( 2-7 

q 10 P 10 

Here X is determined from the water-vapor sounding. In connection with 

the present study, the Northern Hemisphere soundings tabulated by London 

2 

(1957) have been tested to obtain a best-fit technique between In q and 
In p so that an equation of form 




might be obtained. The values of q(p) corresponding to O = 0.1, .3, .5, 

.7, .9 were used to solve for the profile-parameter X, which varied with 
latitude and season. The correlation coefficient which related In q to 
In p was determined with such high precision that use of (2-8) for extra- 
polation to a = 0.1 was considered feasible when q(p) was available only 
at the levels p = 0.3, .5, .7, .9. 

With q ( 1 ) known, the scaled absorber-mass U n modeled after U Q of 

2 o , 10 

Eq. 2-3 was formulated for use in the model. In some cases values of q^ 
will also be missing but such values are subject to retrieval from Eq. 2-8 
using the mixing ratio values from O ^ 0.5. 

The procedure of water-vapor extrapolation to at least O ^ .3 should 
be considered a part of the model presented here with the relevant analysis 

2 

By LT. R. J. Plante in his 1972-73 Master's thesis. 
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routine attached as a subprogram. The values of are quite small, but 
the transmissivities in the terrestrial spectrum are appreciable in any 
case. 

2(c) Criterion for the existence of clouds at gridpoints. 

In both the solar and terrestrial radiative transfer problems, it is 
essential to know the fractional cloud-cover CL at each gridpoint. This 
model, like that of FNWC, permits only a single cloud layer to form. The 
criterion for the amount CL of cloud-cover in existence at any gridpoint 
is given by Smagorinsky * s (1960) algorithm (applicable at approximately 
700 mb): 

CL = 2.0 ( J ) c= ? - 0.7 (2-9) 

s 

Thus at every gridpoint there is a cloud-covered area CL of fractional 
amount 0 < CL ^ 1.0, and a clear-area fraction 0 5* 1 - CL ^ 1.0. 

When the fractional cloud-cover CL exists, the cloud is assumed to 
extend everywhere within the grid area square between a = .8 to O = 0.4. 

The occurrence of CL > 0 of cloud amount profoundly affects the treat- 
ment of both the solar (Sec. 3 ) and terrestrial (Sec. 4 ) radiative 

transfer problems. 

An improved cloud model permitting the introduction of a two- layer 
cloud system both of which fulfill Smagorinsky-type requirements is en- 
visaged for future model development. Such improvements would follow a 
procedure somewhat similar to that of Oliger et al. (1970). This refined 
radiative cloud treatment would permit CL^ of intermediate cloud between 
CJ=.8 and G=.6, and CL^ of high cloud between G=.4 and G=.2. In the model 
of Oliger et al., the high-cloud layer is ascribed a solar transmissivity 
(or alternatively, a reflectivity) intermediate between moist air and a 
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mid-level cloud deck. Moreover, the distribution of the upper level 
cloud elements is assumed to overlap the lower-level ones in a random 
fashion. 

2(d) Integrated absorber-mass formulas 

If the result of (2-2) is summed from p(2k) = (.2k)p Q to p = p^ 
using straightforward integration, one obtains 

C (2k) = 3.14 X 10" 4 — \ i 2 -)' 72 

®>o/> P o p o 

(♦2k)p o 



C(2k) " 3,1 1.72 10 " H[1 ‘ C^ 1 ) 1 * 72 ] (k=0, .... ,5) 



( 2 - 10 ) 



where H = 7.995 x 10 cm is the height of the homogeneous atmosphere. 
Since the surface pressure TT may differ considerably from p Q , it is 
convenient to adapt Eq. 2-10 to the a - levels of Fig. 1. This is 
easily done by setting the ratio 



p ( 2k) = n a 



(2k) 



Thus C(2k) becomes 



.-4. 



C (2k) = - - - 1 - 1 X y^ 0 H (^ ) 1 * 72 [1 - a 1,72 ] 

p o 



( 2 - 11 ) 



where O = 0.0, 0.2,..., 0.8, 1.0 spans the 5-layers shown in Fig. 1. 
For water vapor, a similar integration of Eq. (2-1) over the layer 
(10,2k) gives the reduced absorber mass in this layer: 



xr i o P n=k£l - _ 

i t / oi- \ /TT x 1 * 72 / O ^ v"i " r ~ 1.72 

u(2k) ■ <? > <TT2i> E , «2n-l °2n * a 2 n-2 

o n=5 



1.72 



(2-12) 



In (2-12) ^20-1 tke var i- a ble ffl ean mixing ratio of water vapor in 

1 72 

the layer (2n-2,2n). The parameter (TT/p^) * is constant at a given 
gridpoint. 
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The computation of C(2k) and U(2k) then depends upon the values of 
1*72 

^2k * at ^ eve ^ ^k at eac ^ gridpoint. The following tabular values 
1.72 

of a * are useful in these computations 



Level 2k 


0 


2 


4 


6 


8 


10 


CM 

D 


0.0 


.2 


.4 


.6 


.8 


1.0 


„ 1.72 
2k 


0.0 


.06277 


.20680 


.41535 


.68126 


1.0 


With U(2k) and C(2k) 


given by 


(2-12, 


2-11), 


these mass 


parameters increase 



from the earth 1 s surface upward. 
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3. The solar-radiative heating package. 

(a) Partition of the solar spectrum at earth. 

Following Joseph (1966), a solar constant of 2.00 ly min ^ is assumed, 
but is subject to 4$ attenuation by oxygen and ozone absorption above the 
tropopause. Thus, an effective solar constant 

S = 1.92 ly min 

has been assumed entering the level O = 0.2. 

The effective insolation on any given Julian date D at level 2, Fig. 1, 
is then obtained from 



F 



r -2 

= S (— ) cos 
r 



(3-1) 



m 

Here r/r is listed as a function of the date D in Smithsonian Table 169 
m 

(List, 1963). In addition, the zenith angle 0 is available for any given 
time step (GMT), at latitude Cp and longitude X, using a subroutine of the 
FNWC radiation package (Kesel and Winninghoff, 1972, Appendix B). 

Joseph (1966) has devised a computationally simple partition of the 
solar insolation into the part 



F a = 0.349 F. 

A 

subject to water-vapor absorption but not to Rayleigh scattering in the 
atmosphere. The remaining insolation of the partition, denoted as 

F s = 0.651 F 

is subject only to Rayleigh scattering, but not to absorption in a clear , 
moist atmosphere. The presence of cloud-decks introduce cloud-reflecti- 
vities, however, in both the F and F solar-energy partitions. In any 

A o 

case, the F^-energy comprises all wavelengths X < 0.9 (J. m since here water- 
vapor absorption is negligible, while the F^-category comprises solar 
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wavelengths X ^ 0.9 |a m where absorption by water vapor and CO^ bands 
in the near infrared is the dominant attenuation process acting on the 
solar beam. In summary, the partition just devised gives rise to the 
parts 



f a 


= 0.349S 


,r ^ 

(— ) COS 

T* 


0 


(3-2) 






m 






F s 


= 0.651S 


,r . "2 
(— ) cos 
r 

m 


0 


(3-3) 



The component (3-2) contributes to atmospheric heating in any vertical 
column whenever cos 0 > 0. The component (3-3) contributes only to earth- 
surface insolation, after scattering and reflective attenuations in the 
column above earth are considered (and when cos 0 > 0) . 

In some climatological and general circulation investigations, e.g. 
Manabe and Strickler (1964), it has been found useful to compute the 
time-mean cos 0 for use in (3-2, 3-3) for the particular date under con- 
sideration. The following table, after Manabe and Moller (1961) lists 
3 Q 

daily mean values of cos 0 as a function of latitude and time of the year. 
In addition. Table 2 lists the duration of the day in hours for these 
same arguments. 



TABLE 2. The seasonal and latitudinal distributions of the 
length of the daytime given in parts of 24 hours and those of 
the weighted mean values of cos 0. 



°Lat . 


Fractional length daytime 


Weighted mean 


cos 0 






Apr. 


July 


Oct. 


Jan. 


Apr. 


July 


Oct. 


Jan. 


5 


.508 


.517 


.500 


.496 


.625 


.587 


.614 


.591 


15 


.521 


.537 


.492 


.471 


.618 


.601 


.579 


.549 


25 


.533 


.562 


.483 


.450 


.599 


.593 


.524 


.474 


35 


.546 


.596 


.471 


.421 


.558 


.567 


.458 


.393 


45 


.562 


.637 


.454 


.362 


.501 


.521 


.379 


.317 


55 


.596 


.708 


.437 


.321 


.423 


.453 


.282 


.203 


65 


.629 


.837 


.404 


.208 


.345 


.369 


.176 


.106 


75 


.750 


1.000 


.329 


— 


.241 


.311 


.071 


— 


85 


1.000 


1.000 


— 


— 


.168 


.318 


— 


— 



sec 0 is very nearly equal to the reciprocal of cos 0. 
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3(b) Disposition of Fg-insolation 
A. Clear-sky F^-insolation 

Joseph (1966), found that the Rayleigh-scattering reflection values 

at sea level (after Coulson, 1959) could be fitted by least squares. 

The result was a best-fit regression of form 

a (TT, 9) = .085 + .25074 log (- sec 0) (3-4) 

S P o 

where 7T is the surface pressure and p Q = 1013.25 mb. The value of is 
constrained never to exceed an upper limit of 0.99 even for large sec 6. 
3(b)B. Earth-surface albedo 

Another reflective-type parameter to be considered is that of surface 

albedo 0t , assumed constant over bare land as cc = 0.14 (except over 
g g 

deserts'), but variable with respect to cos 0 over the ocean: 

a g (SEA) = MAX | . 06 , .06 + .54(.7 - cos 0)j (3-5) 

The value of & in the Sahara-Arabian desert is also considered a function 
8 

of cos 0 through 

a (DES) = .8 - .6 cos 0 (3-6) 

8 



The gridpoints which lie in these deserts are easily located in the FNWC 
model as corresponding to I ^ 45 and latitude Cp <S arc sin .57. Also 
within the northern snowline (determined from monthly climatological 
surface temperatures T^ ^ - 3°C), is obtained from the empirical 
result (Posey and Clapp, 1964): 



f ^ C LAT " 

a (SNOW) = 0.4 {1.0 + 

o ^ 



< C LAT- 45 -> 2 + < C LAT - 5 -> ' 



(3-7) 
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Here is the number of latitude degrees poleward of the mean northern 

snowline. A similar expression applies for the south-polar ice-snowline. 

These values of (X have been drawn from either the Kesel-Winninghof f (1972) 
8 

or Gates et al (1971) documentations, respectively. 

At any clear-sky gridpoint, the Fg-insolation is subject to the Ray- 
leigh sky reflectivity given by (3.4), and the appropriate earth-surface 

reflectivity CL . Considering the possibility of multiple reflections 
S 

between earth and atmosphere, each of which diverts downward the fraction 

CL of the earth-surface reflectance CL , the retained insolation in the 
S g 

earth is I given by 
so J 

I (10) = F (1-a ) [ l + a a + . . . (a a ) n + ...](i - a ) 

SO ss sg sg g 

or 



V 1 - v a -sl 



W 10 ^ l - a a 



(3-8) 



s g 

3(b)C. Cloudy-sky insolation in the Fg-region 
(i) In the cloud layer 

We make use of the Rayleigh reflection coefficient CL of (3-4) applied 

s 

only over air above the cloud-top at p = .4TT. Since the cloud layer is 
composed cf large water droplets, it is subject primarily to Mie-type or 
diffuse back-scatter which gives rise to the cloud-reflectivity parameter 
R^, (approximately, R^ = 0.5), and also to the downward scattering which 
ultimately must emerge through cloud-base at p = 0 o 8TT. 



The two sky reflectivity factors a (- sec 0) due to air-scattering 

P o 

and R due to cloud-reflectance, go on independently with some downward 
reflection following each cloud-to-sky reflectance. Hence, an effective 
air-cloud transmissivity through the cloud-top similar to that implied by 
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(3-8) applies for F -insolation penetrating the cloud-top: 

b 

4TT 

F [1 - a sec 6) ] [1 - R ] 

s s p C 

F _(4)l = (3-9a) 

- % V 

Since the cloud-penetrating insolation F (4)J is not subject to absorp- 

S Li 

tion within the cloud, but only to forward Mie-scattering, the fraction 
emerging as transmitted through the base is likely to have the mean slant- 
path for diffuse radiation, so that F ^ becomes 



F [1 - <X sec 0 )][1 - R ] 

s s p D C 

F r (8) i 2 (3-9b) 

[1 - a R] 
s C 

Here sec 0^ = 1.66 has replaced sec 0 as in most cases of diffuse-sky 

slant paths (Katayama, 1966). Wherever ol^ appears either in the numerator 

or in the denominator, the value of TT in (3-4) has been replaced by .4TT 

and sec 0 by sec 0 = 1.66. The value F of (3-9b) then gives the trans- 
it S Li 

mitted diffuse beam entering the subcloud layer. In this layer, F i 

S Li 

will be subject to further multiple reflections to be discussed below. 
3(b)C. (ii) Partition of F sC (8)J in the subcloud layer. 

The assumption is made here that Rayleigh scattering by air molecules 
occurs in the air-layer between p = 0.8TT to p = 1.0TT. Furthermore, the 
mass of air involved in this layer is 0.2TT, and the angle 0 is still 
subject to the diffuse slant-path condition arc sec 0 = 5/3. As a result, 
the subcloud scatter-reflectivity Ot^ becomes, by analogy with Eq. 3-4: 

a (8-10) = .085 + .25074 log I) (3-10) 

s p 3 
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The diffuse beam F s £,(8)J which first enters the subcloud layer is 
subjected to a sequence of reflections depicted schematically in Fig. 2. 
Consider the beam B-.G-, associated with the first entry of F _ ( 8 ) i at 
point of the subcloud layer (Fig. 2). The fraction Ct^ is turned up- 
ward during the first transit and the fraction (1 - a ) reaches the 

s 

earth at G« • In addition to a , the further fraction a (1 - a ) is also 
1 s g s' 

turned upward from the first entering beam B^G^ before earth-absorption 

occurs. The beam B^G^ therefore has been decomposed into the parts 

(1 - a ) (1 - a ) F n into earth at G n 

s /v g sC 1 

[a + a (1 - a )]F starts along G- -B 0 

s g s sC b l 2 

The latter part is assumed to reach the point where the subsequent 

multiplicative fraction R starts down from B towards G 0 . The inter- 

L> Z Z 

vening subcloud scattering layer reduces the transmission along B^G 2 by 

the factor (1 - OL ); hence at G 0 , we have the downcoming second diffuse 
s Z 

beam insolation contribution 

[a + a (l - a )] • (l - a ) r F (8)1 

s g s s C sC 

with a further decomposition into the parts 

[a + a (1 - a )] (1 - a ) R a F _ — reflected at G 0 
s g s s C g sC 2 

[a + a (1 - a )](1 - a )(1- a ) R F _ absorbed at G 0 

sgs SgCsC Z 

Continuation of this process of successive reflections between ground 

and cloudbase yields the final cumulative absorption at the earth, denoted 

as I sC (10) 

I s c<l°) - (l-a s )(l-a g )F sC+ [a« g a-a s )](l-a s )(l-a g )R c F sC 

+ [a s « g (l-a s )](l-a s ) 3 (l-a g)R 2 a g F sC + 
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Fig. 2. Disposition of diffuse-sky solar insolation F (8)i penetrating into the sub-cloud layer G 
A sequence of multiple reflections applied to the beam in giving contributions to ground-penetration 
and cloud-penetrations (leakage) is schematically indicated. 



and on through higher terms* Finally, upon simplification, one obtains 



I sC < 10 > - F sC . ( 1 -V a -V + 



F . c tCy yi-v3<i-y(i-yK c 

Cl - R c d-a s ) 2 a g ] 



(3-11) 



In the above accounting, certain small atmospheric reflection terms 
from the subcloud layer to earth or to cloud have been ignored as small 



It is to be stressed that I „ ( 10 ) given by Eq. 3-11 is the absorbed 

S 

insolation at earth, resulting from the F^-insolation of Eqs. 3-3 and 3-9b. 



ability depending upon the surface and cloud albedos, although R = 0.5 



absorption occurs. The weighted mean earth-surface absorption then 
becomes 



where I is given in (3-8), and I -(10) in (3-11). 
so sC 

A related subject of some interest concerns the leakage back to space 
of the subcloud insolation (Eq. 3-9b) after upward cloud-transmittance. 
Cloud- leakage . It is of interest to compute the solar-leakage through 
the cloud base subject to the approximation noted in the paragraph just 
below Eq. 3-11. Reference to Fig. 2 indicates that the sum total beam- 
transmittances at the cloud-base involve solar energy leaks of sequential 
amounts indicated by the following sum of terms: 



compared to the air-transmission factor (l-<X g ). 



The parameter always arises as a modification of Eq. 3-4 depending 

upon the depth of the scattering layer. 0L and R are subject to vari- 
es u 



is assumed as a first approximation. In the F^-wavelengths no atmospheric 



I (10) = (1 - CL) I + CL I (10) 

S SO SL 



(3-12) 



F 



sC,B 



t = F 



sCl'V V'-V 1 i ( 1 - R c) + E c a g ( 1 - a s )2<1 - V 
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or 



? sC b*" F sc [a s + a e ( 1 - a s )Hi-R c )[ 2 

sC,B sC s g s C 1-R a (1 _ a ^ 

C g s 



] 



(3-13) 



Comparison of (3-13) with (3-11) indicates that the leakage is closer 
to the order of magnitude of the second (smaller) term of (3-11). Choos- 
ing OL m 0.05, and average values Co = 0.14, R = 0.5, indicates the leakage 

S g C 

to be about 0.1 F (8)4 or about 3$ of the extraterrestrial insolation of 

S L» 

Eq. 3-1. This leakage represents an additional component to be added to 
the planetary albedo in cloud-covered areas only. Again it is to be 
stressed that none of the F g radiation is absorbed by the atmosphere in 
this model. 

3. (c) Disposition of F -insolation 

A 

(A) Clear areas 

Here absorption of the F^-component of the insolation is considered. 

Use of the Manabe-Mtiller solar absorptivity^ function a(u sec 0) in the 
form 

a = .0946 (U 2 2k sec (3-14) 



will be used extensively in this discussion as well as in that of the 
cloudy-sky case discussed below. The form of Eq. 3-14 indicates that the 
absorption function a is applied over the corrected water vapor path 
between levels 2 and 2k (see Fig. 1) along the zenith-path of slant angle 6. 
Actually (3-14) is applicable to the entire solar energy F. Restricting 
absorption to the ratio 0.349F requires that the multiplicative constant 
4 

The symbol ^denotes the fractional absorptivity function, while A(2,2k) 

denotes the absorbed energy between levels 2 and 2k. 
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of (3-14) be adjusted to (.0946/0.349) = .271. Hence the absorbed insola- 
tion in the layer (2,6) becomes simply 

A(2,6) = .271 (U 2 ^ sec 0) ,3O3 F A (3-15) 

Similarly, the absorbed insolation in the layer (2,10) is 

A(2, 10)= .271(U 2j10 sec 0) * 303 F A (3-16) 

The absorbed insolation in the layer (6,10) is then given by 

A(6, 10)= A(2, 10) - A(2,6) (3-17) 

Finally, the difference formed by subtracting A(2,10) from F. gives the 

A 

direct transmission of F -insolation to the earth's surface. This direct- 

A 

transmission [F A - A(2,10)] is denoted by F A (10) and is partially reflected 

(the fraction Cc is applicable). The remainder (1 — CX )F (10) is absorbed at 
8 8 A 

earth. The F A ~ insolation absorbed at earth is given by 

I a £10)- F A [l-.271(U 2>10 sec 0) * 3 ° 3 ] (1-a ) (3-18) 

In the foregoing derivation, solar absorption was computed only below 
level 2 of Fig. 1, as a 4$ stratospheric reduction of the solar beam was 
previously applied in the layer above level 2. Likewise no accounting of the 
terrestrial spectral absorption in this topmost layer will be made. This 
implies the assumption that the layer .2 is in continual radiative heat 

balance at all gridpoints. 

3(c).B. F -insolation; Cloudy-sky areas 

A 

In the layer (2,4), the Manabe-Moller absorption is employed in a 
manner similar to Eq. 3-15 apart from the differing water-vapor mass ^ 
A(2,4)= 0.271(U 2 4 sec 0)* J (3-19) 

The F -insolation at level 4 is therefore 
A 

F A ( A)= F a [1 - .271(U 2>4 sec 0)* JUJ ] (3-20) 
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Of the energy F A (4), the fraction 1 - R c is transmitted through the 

cloud top which interfaces at level 4. Hence the downward transmitted 

F -insolation through the cloud top is 
A 

F A (A) * = F A (1 " R c )f 1 -- 271 (U 2>4 sec Q) • 303 j (3-21) 

Within the cloud layer, the water-vapor path in layer (4,6) is subject 

to the mean diffuse-path augmentation factor 5/3 (after Katayam a, 1966). 

Thus the energy-absorption A(4,6) is given by 

5 U r 303 

A(4,6) = F A (4)i *.271[|(U 4>6 + (3-22) 

In (3-22), the water-vapor parameter U is the water ~vapor mass equivalent 

of the cloud-droplet absorption. The determination of the parameter U c 

which Is spread vertically over the cloud-layer (4,8) is modeled after 

estimates of cloud-droplet absorptivity a [given by Korb et al.(1956)] 

OLD 

and makes use of the modelling assumption 

a CLD - - 271 { [ f (U 4,8 + V>' 3 ° 3 - [ ! U 4,8 r3<>3 } < 3 ‘ 23 > 

The method of computation of U /U. 0 is detailed in Appendix C as a function 

L 4 , o 

of g for the latitude bands 0-40, 40-65, 65-90. 

With U^, known at each gridpoint, the following three computations of 



insolation are made 






A(4,6) = F A (4)i 


{• 2711 X > 6 + r»' 303 } 


(3-24) 


A(6,10)= F A (4)i 


{• 27 i[f(u 4>10+ u c )]- 303 -[|(u 4i6+ ^)]- 303 } 


(3-25) 


r gc 4 =f a< 4 >‘ 


{l-0.271[f(U 4>10+ U c )l- 303 } 


(3-26) 
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These three computations represent respectively (i) the absorbed insola- 
tion in the layer (4,6), (ii) the absorption in layer (6,10), and (iii) 
the transmission to the earth's surface. The last of these is subject 
to the multiplicative earth-absorption factor P (of form similar to 
Eq. 3-8) 

P = (1 - <0/(1 - a R ) (3-27) 

& 6 ^ 

This factor P when multiplied by I J of (3-26) gives the F^-penetration 
of the earth's surface under cloudy skies. 

In summary, the gridpoint absorption of F -energy in earth now becomes 

A 

the cloud -weight ed sum 



I A (10) = (CL) (P of 3-27) (I c of 3-26)+(l - CL)(I^ 0 of 3-18) (3-28) 

Similar cloud-weighted sums of absorption of F -energy should be formed 

A 

for the layers (6,10) and (2,6). In the latter connection, it should be 
recalled that the F^ cloudy-sky result in layer (2,6) is made up of the 
two contributions (3-19), (3-22), the former of which has no cloud-reflec- 
tivity factor (1 - R ) since it pertains to A(2,4) only. 

3(d) Summary of solar-heating effects 

The results of Sections 3b, 3c are now combinedl to get a resultant 
insolational-absorpt ion at earth. This is done by summing the two contri- 
butions I (10) of Eq. 3-12 and I (10) of (3-28). 

b A 

Within the atmospheric layers (2,6), (6,10), only the F^ wavelengths 

contribute to absorption and therefore to temperature-change through 

2 



ST,« A(2,6) cals/cm /hr 

SF (2 ’ 6) = , 1/1000 : 



(3-29) 



.4n(— )[c p ] 
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(3-30) 



dT, A 1fn _ A(6 , 10) cal /cm /hr 

St <6 - 10) - , 4n( 1000 

8 P 



The first of these temperature-changes involves the cloud-weighted com- 
posite absorption rate formed from A(2,6) of (3-15) and from A(2,4), 
A(4-6) of Eqs (3-19,22). Similar weighted means are needed for A(6,10). 
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4. Formulation of net cooling rates by terrestrial radiation 
(a) Clear-sky areas 

It should be recalled that Table 1 contains the format of a "radiative" 
sounding having input data defined at the levels listed in Table 3. 



Table 3. Levels and data used in the "radiative a 
soundings . 



Levels 


10 




8 


6 


4 


2 


1 


0 


Temperature 


• 

T 10 = 


T s 


T 8 


T 6 


T 4 


T 2 


T 1 


T o 


W. V.mass 


U 10 = 


0 


U 8 


U 6 


U 4 


U 2 


U 1 


U 0 


CO 2 mass 


C 10 = 


0 


C 8 


C 6 


C 4 


C 2 


C 1 


c o 



The value of water-vapor reduced absorber-mass - 5,....,0) was defined 



in Section 2(d), as was also the corresponding value of C These values 
have been incremented upward relative to the ground as indicated by Eqs. 
2-12, 2-11, respectively. In these equations the notation U(2k) has been 
treated interchangeably with the more compact notation used here, and 
in the remainder of this section. 

Consider the radiative temperature change in air. By the first law of 
thermodynamics, this is given by 

* V 

- F, J 

(4-1) 



ST < F 6 - 



St o. 4 tt 2222 c 
g P 



for the layer (6,10). 
given by 



For the layer (2,6), 



Bt 

ST = • 




0.417 



1000 

g 



c 

p 



the corresponding formula is 



(4-2) 



23 



The numerators of (4-1), (4-2) which are contained within parentheses 
are understood to be flux-divergences from the atmospheric columns 
(6,10) and (2,6) over the gridpoint having the indicated "radiative" 
sounding. 

Virtually all of radiation packages in the presently used large- 
scale circulation models in the United States [cf., Manabe and Strickler 
(1964), Oliger et al. (1970), Arakawa et al. (1969)] make use of Yama- 
moto^ (1952) water vapor flux-emissivities as a basis for their radia- 
tive flux calculations. These absorptivities £(U,T) as defined below 
are virtually independent of temperature in the normal range of tempera- 
tures (300K ^ T ^ 220K) , but at temperatures T < 220°K the wave-averaged 
f lux-emissivity function defined in (4-3) after Yamamoto (1952) and 
others, has the form 

€(U,T < 220) = 2 Hi 

V dB v 

j (V » T) dv 

o 

It is seen that the numerator of (4-3) has its temperature-dependent terms 

3 

fixed at 220K, while the denominator equals 4St*T even when T < 220K. 

Thus (4-3) introduces a temperature-dependent f lux-emissivity for tempera- 
tures T < 220K, reflecting the effect of increased absorptivity by the 
rotational band of atmospheric water-vapor which takes effect in view of 
Wien's law at these low temperatures. In the discussion which follows, 
the absorptivity will be restricted to be due to water-vapor alone, until 
the subject of is specifically introduced. Although analytic functions 



dB v 

[ 1- T (L U)] —T (T=220K)dV 



(4-3) 
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have been determined for £(U,T) for both water vapor and ^ use f u l 

to examine schematic versions of the Yamamoto method of net flux calcula- 
tions through key levels of the sounding (Table 3). Schematic net flux 
depictions are shown in Fig. 3 for the levels 10., 6., 2. in cases (A), 

(B), (C), respectively. 

4a. (A) Net flux at level 10. 

Sasamori (1968) has empirically recomputed e(U) after Yamamoto’s (1952) 
values for the temperature-independent part of atmospheric radiational 
soundings, schematic versions of which are shown in Fig. 3. The tempera- 
ture-independent emissivity G(U), after Sasamori (1968) is 

e(U) = .240 log 10 (U + .010) + .622 (4-4) 

in the temperature- independent region T ^ 220K. For T < 220K, Sasamori 
gives for the temperature-dependent emissivity 

£(U,T) = (8.34T- 353lO8 10 U - 44 ) D -°34551og 10 D-.705 (4-5) 

Here U represents the final U-value relative to the reference level 
considered to be level 10, so that 



i(U,T) = (8.34T ,3531 ° 8 10 U o"* 4A ) U q ~ * 03 ^ 551og io U o” * 75 



(4-6) 



Corresponding to this value of U q , a temperature-dependent 6-isopleth 
converges into the lower left corner of each panel of Fig. 3, i.e., 6 “* 1.0 
as B(T=0)= 0, where B = St*T 4 is the Stefan-Boltzmann flux. 

The trapezoidal summation-rule gives a representation of the hatched 
(downward flux) area of Fig. 3A as 



Fj(10) = .5 |e(8-10)(B 10 -Bg) + [ e(8-10)+ e(6-10)](B g - Bg) 

+ [6(6-10)4- 6(4-10)] (B,-B,)4-[6(4-10)4- 6(2-10)](B -B ) 

° 4 Bdj) 

+ [e(l-10)+ e(2-10)](B 2 -B 1 )+ 2^ 6(0-10, T)dBj- (4-7) 

B=0 
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where e in the final term is the function of U,T listed in (4-5). The 

4 

notation 6(8-10) denotes e (U Q -U ir .) etc, while B in = St*T.. n , 

o 1U 1U 

St = 5.67332 x 10 ^cal cm ^ hr * °K ^ and so on for all of the other black- 
body fluxes (Bg, Bg, B^, B^, B^, B^), distributed along the sounding of 
Table 3. 

Note that the final integral \ £ dB of (4-6) may be evaluated by the 
integration 

T=T 

\ e dB = 8.34 U Q ^■• 03 ^ 551 °gio U 0'* 705 ^ T* 3531 ° 8 10 U 0"* 44 4St*T 3 dT 

0 



= St*T 4 1 {8.34U 0 - 034551 ° gU 0- 705 . 



4t . 3531ogU Q -.44^ 
.3531og 10 U 0 +356J' 



It finally follows that 

. B i 

) ' £ dB = “V * £(U 0 >T l ) * .3531og 4 oV 3. 
o 



56 



(4-8) 



with £(Uq,T-^) defined in (4-5) 
* 

Hence F^q is given by 



F 10 - B 10 - V 10 > 



(4-9) 



with the latter term given by Eq. 4-6, 4-7. 

B. Net flux at level 6 

In this case, the desired net flux is represented by B(T^q) minus the 
slant-hatched area above the sounding curves of Fig. 3B. Thus the result 
may be modeled after the results 4-6, and 4-7 

-.5 {£(6-8)(B g -B 6 )+[£(6-8)+ £(6-10) ] (B 1Q - B g ) 

+ £(4-6)(B 6 -B 4 )+[£(4-6)+ £(2-6)](B 4 - B 2 ) 



* 

F = B 
6 10 



+[£(2-6)+ £(1-6) ] (B 2 - B x ) 



+ (St)T x * e < U 0 > 6’ T l ) * .353 logU Q 6 + 3.56} 



(4-10) 
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Here €(U,T) is defined in (4-5) while £( ~ |) *- s defined i n (4-4). 

C. Net flux at level 2 



Fig. 3C is relevant for this computation. The net flux at level 2 
may therefore be written as follows by analog with Eqs. 4-8, 4-6, 4-5, 
4-4: 

F 2 *= St*T 10 4 - .5 |e(2-4)(B 4 - B 2 ) + [€(2-4) + 6(2-6)] (B fi - B 4 ) 

+ [€(2-6) + 6(2-8) ](B g - B fi ) + [€(2-8) + 6(2-10) ] (B l0 - B g ) 

8.0 



+ CdOZHBj-Bj) + St*^ 4 * eC0 0>2 T 1 )* .353 3 . 56 J 



(4-11) 



where B^ = St*T. 

k k k 

With F 0 , F, , F in formulated in terms of radiative-sounding parameters 
Z b lu 

by Eqs. 4-7 and 4-9, 4-10 and 4-11 respectively, layer-mean temperature- change 

rates in °K(hr) ^may now be calculated. Before this is done, however, it is 

* k k 

convenient to introduce the modifications in , F^ , F^ resulting from the 
inclusion of CO^ flux-emission. 

k 

D. Inclusion of CO absorber-mass in F calculations 

r U N 

The "area” integrals \ C(U,T) dB depicted in Fig. 3 have been computed 

U=0 

by integrating along a two-path sounding curve away from the reference level 
(Fig. 3A,B,C). When absorption is due both to water vapor and C0^ 9 the 
product-law of transmissivities must be employed. In doing this the trans- 
missivity of the two absorbers acting jointly is expressed as 

T f (U,C) = T f (U) T f (C) 



Then the flux-integral becomes 

U , C 

AT* 



_V C N 



l 



N* N f 



F* - ^ [l-T p (U)T F (C)]dB = ^ |[1 -T f (U)] + T p (U) e p (C)} 

U,c=0 U=0,C=0 "w.v.flux" "CO " 



dB 



(14-12) 
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The first term on the right is the water-vapor flux computed using for 



example Eq. 4-7 and ignoring the emission by CO^. The second-term, 
denoted "CC^", represents the effective emission due to CO^ not already 
counted in the water-vapor flux computation* The added emission is 
expressed schematically on the Yamamoto chart by an increase of emissivity 
at all points of the radiative soundings (except at the reference level 
points). Thus at any temperature T, the required joint-emissivity is 
obtained by replacing that due to water-vapor already computed, €^(U) , 
by 



€ 

WC 



= V u > + A€ wc 



(U,C). 



The last term includes the 15 fim band emitted flux of CO^ which has not 
already been counted in emission by the rotational band of water vapor. 
Sasamori (1968) has computed this 15 band "overlap" region emissivity, 
and developed the empirical correction shown in (4-13, 4-14) below: 



i€ W,C <U 2k> C 2k> - C (» a ) * S c (C 2k ) 

Ae wc (2k) = .07262 jl-. 62556 (U 2k +. 286) * 26 j {log C 2k + 1 



. 064 } 



(4-13) 

(4-14) 



Thus in considering the flux-modification for the inclusion of CO^, 
one simply computes the added emissivity-values of (4-14) at each 

sounding point. In the summation Eq. 4-7, for example, where e(8-10) 
appeared, this term is simply replaced by the function 

V U 8 - U 10> + A€ WC [( V V’ <C 8- V 1 

where the function ^ e wc term uses the simultaneous values Ug - U 10 and 

Cq — C- stored from the radiative sounding. The 15 (Jm effective emissivity 
8 10 

Ae wc of Eq 4-14 is considered to be temperature-independent (Sasamori, 1968) 
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in the range T ^ T 1 ^220K, in which range the water-vapor emissivity € 

1 w 

was temperature-averaged. Hence in Eq, 4-7, the final term of F*(10) 
becomes 



St * T 1 te w (U 0 ,T l ) ( .3531og 10 U 0 +3.56 



) + e wc of (4-14) J 



(4-15) 



where Ae _ in the last term of the bracket in (4-15) involves substitu- 
WC 

tion of Uq, Cq for U,C in (4-14). 

* * * 

With the revised formulations for F 2 > Fg , F^q it is appropriate 
to compute the flux divergences in the layers (2,6), (6,10) and finally 

F io* - B 10 - F W c< 10 >' (4 - 16) 



which gives the terrestrial radiative loss at the surface with clear skies 
above. 

In addition, only a simple programming alteration is required to 

★ * 

include the CO^ modifications to recover the flux divergences - F^q 
jc * 

and F 0 - F, necessary for the computation of the layer temperature- 
Z o 

change rates, Eqs. 4-1, 4-2. 

* 

4(b) Modification of net-flux F with cloudy skies. 

As before, the cloud-cover amount CL is given by Eq. 2-9. Section 4(a) 

has dealt with the clear-sky fraction (1-CL) over gridpoint (i,j), and 

the remaining fraction CL therefore represents an overcast (extending 

vertically from O = .8 to O = .4) over the remainder of the grid area. 

★ ★ ★ 

As before, we are interested in F^q , Fg , (but with overcast) over 

the gridpoint. The relevant schematic flux diagrams, after Yamamoto 
(1952), are depicted in Fig. 4. In Fig. 4, all emissivities are based upon 
the function as computed using Eqs. 4-4 and 4-14 in proper combination 
of absorber masses relative to the reference level. 
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Fig. 4. Schematic area-depictions of net flux at reference levels 10 

and 2, when an overcast cloud extends between levels 8 and 4 in the 

vertical. ^The horizontally-hatched area is ; the slant-hatched 

area is F . The net flux F * is zero, since it lies within the black- 

body cloud. F * is numerically equal to F in . 

o 1U 
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4(b)A. Net flux F^ Q at earth, with overcast clouds. 

The result for this case is based upon Fig. 4 and is given by the 
horizontally-hatched area: 

F 10 = (B l(f B 8 )[1 ■ * 5 e wC (U 8-10 ,C 8-10 )J (4 ” 17) 



This is a cooling-rate at the earth-surface. 

At the same time the net flux at level 8 (= the overcast-base) is an 

identically equal magnitude, so that the layer 8-10 is in radiative equilibrium. 

However, the cloud base receives net flux warming exactly equal to 
* 



f 8 ® - F 10 > a warming 



(4-18) 



4(b)B. Net flux at level 6. 

Since a layer of approximately 50 meters of average cloud-drop content 
may be regarded as a black body (Brunt, 1939) at the mean temperature of 
the cloud layer, it follows that level 6 is characterized by 

F 6 * = 0 (4-19) 



4(b)C. Net flux at level 2. 

The design of the result in this case is suggested by Fig. 4, 

*• 

F^ becomes 

F 2* " V • 5 i e wC (2 - 4)< V B 2 )+ C «C (1 ' 2)(B 2' B 1 ) + 
4.oe(u T,) 

2B 1 [353 log 10 U 0 2 +3.56 + Ae wc (u 0,2 C 0,2 ) ° f (4 " 14)1 } 



Thus 



(4-20) 



after application of Eqs. 4-11, 4-12 to the truncated area of Fig. 4. 

This net flux is a much smaller value than that of Fig. 3C because of the 
cold black body interposed at level The reduced net flux due to clouds 
is an effect often referred to as "cloud-contamination" in problems of 
remote sensing of the atmosphere. 
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4(c) Terrestrial radiative cooling rates in key atmospheric layers 

(i) Upper layer (2,6 ) 

The radiative flux divergence in the cloud-covered case is 

(F 2 * of 4-20) - F 6 *(= 0) (4-21) 



This result, appropriately cloud -weighted with the clear sky result 
★ * 

F 0 - F, of Section 4(a), then gives the terrestrial radiative cooling 
c o 

ic 

in the layer (2,6). The cloud-weighted - F^ may then be used in 
connection with Eq. 4-2 to yield a temperature-change rate due to 
terrestrial radiation. The latter may then be combined with the insola- 
tional absorption of 3-29 to obtain an overall radiative change in the 
layer. 

(ii) Lower layer (6,10) 

The corresponding result for this cloud-covered layer is 

( P 6 * = 0) - Fg* = Fg* = - F 10 * (4-22) 



where a negative result indicates a net flux-convergence, or a warming 
rate. The result of (4-22) covers the layer (6,8). To this may be added 



the zero-change in layer (8,10), so that 

F 6 " F 10 = ~ F 10 ° f ^~ 17 ^ 

(iii) The interface level a = 1.0 

The terrestrial radiative transfer across the earth-surface is F 
which must now be cloud-weighted in the usual manner, i.e. 



(4-23) 



10 



F 10 *(i,j) = CL [F 1Q * of (4-17)] + (1-CL) [F 1() * of (4-16)] 



(4-24) 
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5. Maintenance of the heat budget 
5(a) In the atmosphere 

At each of the odd levels in Fig. 1, the diabatic heat-sources of the 
radiative type discussed in Sections 3 and 4 are included in the thermo- 
dynamic equation 

St (nc P V + v • <vA> + S3 'ViV - A 

. 

Here is the diabatic heating (cooling) rate per gram at level k includ- 
ing that due to radiation. should also include the latent heat released 

at level k at time-step t. The subject of latent- heat release and more 
generally of the atmospheric water vapor balance is not included in this 
paper. 

5(b) At the earth* s surface 

The earth* s surface is considered to be a heat reservoir of zero heat 
capacity over land, and of infinite heat capacity (a perfect insulator) 
over the ocean. The case of a solid ice sheet may be considered to be a 
predictable compromise between the land versus ocean cases. These earth- 
surface vertical heat-flow requirements lead to the following surface 
boundary conditions, essentially after Langlois and Kwok (1969) henceforth 
abbreviated to LK. 

5(b) A. Over land 

The external radiative surplus (deficit) at level 10 is 

* 

I - F 
10 10 

With the land surface considered to be of zero heat capacity, this surplus 
must be diffused upward in the form of sensible-heat conduction T and 
latent heat of evaporation E at the earth surface, across land: 

X 10 " F 10* = r +E=r[l+i] (5-1) 
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In (5-1), r - F/E is the so-called Bowen ratio, given empirically by Sellers 
as 



r = 9.6 (sin Cp) 2 - 7.93 (sin cp) + 2.0 (5-2) 

sin Cp ^ 0*4 for all Cp 

In this treatment, a direct formulation of E over land is not achieved; 
rather, only an indirect computation making use of E = T/ r is employed. 

The determination of the sensible-heat flux T is based upon its large- 
scale parameterization [(after Deardorff 1966), LK(1969)], and a simultaneous 
micrometeorological parameterization. Thus LK give the large-scale para- 
meterization as 

r "PlO C p K Az 



r = 




V 9 



=■5 



(5-3) 



and the micrometeorological equivalent as 

r * P 10 C p C D V S (T g - V 



(5-4) 



In (5-4), V g is related to the surface geostrophic wind through V g = *8^1 
+2.2 in mps. In (5-3) K , a and 7 C are turbulent-transfer coefficients 
whose values have recently been documented by Kaitala (1972). The value 



of 0g is related to the estimated T^q = T^ by 

a m T /JL» r 286 

D s lonooo' 



(5-5) 



where T - T is the actual temperature lapse (compatible with Eq. 5-3) within 
g x 

the constant flux boundary layer. 

A(i) Overland boundary adjustment of surface temperature 

The combination of Eqs 5-3, 5-4 leads to the solution for T x at the 
top of the boundary layer: 
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(5-6) 



( V Z 10 } [C P V S T g - Ky C ] + K 9 9 

T x- [K( 10 ffi ) . 286 +CDVg(v Zio)] 

In (5-6) 0g is related to by a formula of form (5-5). If one starts over 

land with T = T„, a solution for T emerges and a value of T is available 
g S x 

from (5-4). Application of the Bowen-ratio technique then provides a value 
of E. With first iterates of T and E available for use in (5-1), an adjust- 
ment in surface temperature is achieved as follows 



F io*<y ■ *10 - r(1 + b 



(5-7) 



F 10 B 10 " F d F 10 <V + ^ T g " V 

u 



(5-8) 



If one makes use of Eqs 5-7, 5-8 for correction of T in 5-6, a revised 

r of (5-5) results and sequentially there is a second iterate T j from 5-7, 

5-8. This process must converge rapidly in view of the large coefficient 
3 3 

4 Tg or 4 T^ , depending upon the step in the iteration procedure. When 

the several iterative steps are complete, a heat balance exists at the 

earth* s (land) surface, which has realistically determined the surface tem- 

: k 

perature. Induced changes in F^ and have resulted from the recalcula- 
tion of the earth-surface temperature to be T^(instead of Tg). These become, 



respectively, in the case of clear-skies 

= [B(T )- B(T g )J [1 - e wc (6-10)] 

AF 2 * = [B(T gn ) - B(T S )] [1 - € wc (2-10)] 

where n is the iterate-step and is given by Eqs, (4-4, 4-14) 



(5-9) 

(5-10) 



With overcast skies only the minor correction below applicable in the layer 
(8,10) occurs 

AF*(8,10) = [B(T ) - B (Tg) ] [1 - e wc (Ug >1() )] (5-11) 
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A(ii) Over the ice-free ocean. 

Here it is assumed that Eqs 5-5, 5-6 still apply, but that T is known 

8 

from the sea-surface analysis, which is held constant during the prediction 
process. In this case, a unique value of and therefore of T results. 

It is no longer necessary to make the Bowen ratio assumption for the surface- 
layer evaporative heat-transfer rate. Instead the analog of (5-3) is formed 



E = L Pio 



K* 



* 0 9 
1 + a (■ 



Z 9 _Z 10 




and is coupled with the analog of (5-4) 

E * L PioVs (< isg- V 



(5-12) 



(5-13) 



Eqs. 5-12, 5-13 yield the result 

■S + <V Z 10> C D V S q Sg 
q * " K + (2,- Z 10 > C D V s 



(5-14) 



with 



.622 e (AE - BE/ V 



S§ n - e (AE-BE/T g ) 



(5-15) 



AE = 21.656; BE = 5418. The last relationship is the Clausius-Clapeyron 
equation. 



In this case (5-1) may be rewritten 



* 



F 10 ■ ho • <E + r > 



(5-16) 



If of (5-16) differs from that computed using Eq. (4-16), the differ 

ence I - F^* should be used to increment both E and T in the relative 
proportion of the magnitudes of the latter pair. 
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A(iii) Over ice-covered ocean 

In this case, the surface-to-air boundary condition may be stated as 



r = I io- F io - b < T w-y 

-2 -1 -1 

where b » ice-conduction coefficient (697.8 erg cm sec deg ) 



(5-17) 



= 271. 2K is the freezing point of sea water 



T = 
g 



the desired temperature at the ice-air interface 

The evaporative-heat transfer E has been omitted from 5-17. However, Eqs. 

5-3, 5-4 are still valid and it is still necessary to obtain in order to 

evaluate T by 5-4. Accepting T = T as the first iterate, it is possible 

g ^ 

* 

to adjust F 10 (T ) of (5-17) as was done in Eq. 5-8. If only a single 
iteration is required, an analytic formula for results [Kaitala (1972); 



LK(1969)], namely: 



CD + 



pe c 



T = 
x 



Z 9" Z 10 



- 07 , 



Z -Z„ \9TT 
9 S 



( 1000 ) .286 | 



p 10 C p C D V S (1 < ) 



(5-18) 



where 



D - ho + bT w + 4St * T S 4 - F 10* (T S) 



- »10 C P 



K* 



u ♦ 



V 2 S 



£ = (^/(Cp + b + 4St*T s J ) 

°r = Pio c P Vs 



(5-19) 
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With known in terms of T^, one may iterate successively 

values of the latter using 5-8 and obtain convergent values of 

well as of the sensible heat transport T. With T replaced by 

S 

realistic value, T , it is relevant again to replace 

S 

* * 



following the method 



F 6 by 


6 g 


* 


* 


F 2 by 


F 2 (T g > 


F io* by 


F 10* (T g ) 


of Eqs. 


5-9, 5-10 



upon the 

T as 
g 

a more 
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APPENDIX A 



Vertical extrapolation of the mixing ratio profile. 



When London's (1957) climatological soundings included mixing ratios 
q(p) to levels of p = 0.3TT, the interpolated values of q(o) (where 10. 0a = 

9., 7., 5., 3.) were subjected to a best-fit against pressure in the loga- 
rithmic form 

In = X in 2- = X In - ( A _l) 

q 5 P 5 .5 



Eq. (A-l) was tested for eight latitude belts on the Northern Hemisphere and 
over four seasonal mean-soundings at each belt. The parameter X is related 
to the correlation coefficient r(Y,X) (where Y = In q/q^ and X * In p/p^), 
which was well in excess of r ^ 0.99 in each sounding. This indicated an 
effective method of extrapolation in the form 



S_ = 
q 5 




(A-2) 



as also suggested by Smith (1966). This relationship permitted extrapolation 
of individual moisture soundings to p = .3TT when q(p) values existed at 
p = .9TT, .7 tt, • 5TT at analysis points. 

When moisture values q(p) exist also at p « • ITT, a relationship of form 
(A-l) was also found to give close verification at o =* O.lTT in terms of a 
best-fit power law 



a = <“> 

q 3 p 3 



(A-3) 



with X determined using the four pairs of (q/q^* p/p-j) values from the 
analysis covering each gridpoint. 



40 



APPENDIX B 



Solar undepleted insolation in spherical coordinates. 

The incoming undepleted insolation is given by 

F = S[sin Cp sin 6 + cos Cp cos 6 cos h] (— )^ (B-l) 

m 

where (r/r ) is listed by Julian date in the Smithsonian Tables (List, 1963) 
m 

number 169. Likewise, the declination 6 is listed in that same source, but 
may also be computed from 

sin 6 = .39785 sin[4. 88578 + .0172D + .03342 sin(.0172D) 

-.001388 cos(.0172D) + .000348 sin(.0344D) -.000028 cos (.0344D) ] (B-2) 
The hour angle h on day D is given in terms of GMT (denoted by Z) by 
means of the relationships 

t = M0D[(Z + 2 ^ 724 ) + 24 » 24 1 

cos h(Z,\) = cos (t- 6 ) - j] (B-3) 

The expression in brackets in (B-l) is cos 6 , where 8 is solar zenith 
angle at latitude Cp, 6 and hour angle h, the latter given by B-3. If 
cos 8 of (B-l) becomes negative, it is set to zero. 
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APPENDIX C 

Determination of the absorber- equivalent water vapor 
mass for droplet-absorption in clouds 

In clouds, the water droplet-absorption has been modeled using the dif- 
fuse-path absorptivity of Manabe-Moller (cf., a^p °f 3-23). In this paper 
the value of a rT , due solely to cloud droplets has been taken after 
theoretical computations of Korb, Michalowski and Moller (1956), and sub- 
sequently used also by Manabe and Strickler (1964). These values are as 
listed below in Table 4. 



Table 4. Absorptivities (a rT ) of cloud droplets 

CLiD 



Cloud-type 


a CLD 


R 

c 


Equivalent Lat. band 


Low 


.035 


.21 


0-39N 


Middle 


.020 


.48 


40-64N 


High 


.005 


.69 


65-90N 



Since only a single cloud-layer has been invoked in the radiative model 
presented here, the larger absorptivities associated with low clouds (lar- 
gest droplets) have been associated with an equivalent latitude band 0-39N; 
while intermediate values of a have been identified with an intermediate 
latitude band, 40-64N. Finally, the minimum value a^^ has been ascribed 

to the latitude band 65-90N. The values a _ of Table 4 are due solely to 

ULD 

absorption by either the liquid-water or the solid-ice content, depending 
upon the cloud type considered. The values a p of Table 4 are then to be 
considered absorption-means for the latitude belts under consideration. 

No specific use has been made of the values of R^, although it seems reason- 
able to ascribe these reflectivity values to the same latitude bands for 
which a CLD -values are proposed in Table 4. 
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The absorptivities of Table 4 apply over the full energy range of the 
solar spectrum. If cloud-water absorptivities are to be modeled as water- 



vapor equivalents occurring solely in the F^-wavelengths , then each value 
of a must be multiplied by (1./.349) on the left side of (3-23). Thus 

C/tiJL) 

one obtains the water vapor equivalent of the cloud-water absorptivity in 



the form 



a* = 2.8653a CLD = -j [|(U 



1‘3^4,8+V 1 ' 303 - 'K, 8 r303 }<- 271) 



(C-1) 



using a Manabe-Moller parameterization similar to (3-23). In (C-1), the 
absorption in clouds has been spread over the full cloud-depth which con- 
tains U, Q , a known absorber mass of water vapor as determined from the 
4,8 

sounding. Eq. C-1 may be solved for U^/U^ g using the appropriate value 
of a = 2.8653a and Table 4. The desired solution assumes the form 

U 4,8 U 4,8 



which becomes 
U 

c 

U 



4,8 



l [1 



.271 ''U 



.6 . 303 j 



.303 



4,8 






(C-2) 



where the exponent 1^0.303 = 3.3003 

A graphical solution of (C-2) for as a function of g is presented 
in Table 5 as a* varies over the low, middle and high cloud (i.e. latitude) 
types of Table 4. 
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U 4,8 


u c (D 


U c (2) 


U c (3) 


0 .Cl 


0.14170 


3.05110 


0.00740 


0.02 


0.18936 


0 .0 7292 


0 . 3 1155 


0.00 


0.26c 06 


0. 12008 


0. 02103 


0. 10 


0 . 4 C 8 5 2 


9.17651 


0 . 0 _> 3 2 9 


0.15 


0 .50529 


0.22673 


C. 04365 


0.20 


0. 53595 


0. 2 o939 


0.05294 


0.30 


0 . t i 7 0 5 


0.344 73 


0.06955 


0.50 


0.9 86 17 


0.47267 


0 . 0 4 B 2 6 


0.80 


1. 2574C 


0.63548 


0 .13520 


1 .00 


1.401 1 9 


0. 73220 


0. 15739 


1. 20 


1.65219 


). 8 22 7 1 


0.17822 


1.40 


1 .31232 


0 . V 0 8 2 


0.198 30 


l . i o 


1. GO 6 48 


0. 98574 


0.21691 


1.80 


2.11306 


1.06793 


0.235 10 


2.00 


2.29408 


1 . 14331 


0. 2 62 o 7 


2. 20 


2.39016 


1.21623 


0.26970 


2. 40 


2.62206 


1 .23659 


0.28626 


2.60 


2.65022 


1. j5504 


0.30240 


2.80 


2 .77499 


1.42297 


0.31816 


3.00 


2.05675 


1.4 3o54 


0. 33357 



Table 5. Water-vapor equivalent U for droplet absorption in 
relation to the pressure-scaled wafer-vapor absorber mass g 
(in the sounding layer (4,8)). Case U (1) corresponds to * 
low-cloud water-drop absorptivity; D (2) to mid-level 

cloud-drop absorptivity; U £ (3) to high-cl8ud absorptivity. 
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